Dynamical structure factor of the anisotropic Heisenberg chain in a transverse field 



m 
o 
o 



Jean-Sebastien Caux 1 , Fabian H.L. Essler 2 and Ute Low 3 
ITFA, U. of Amsterdam, Valckenierstraat 65, 
1018 XE Amsterdam, The Netherlands 
2 Department of Physics, 
Brookhaven National Laboratory, 
Upton, NY 11973-5000, USA 
3 Institut fur Theoretische Physik, Universitat zu Koln, 
Ziilpicher Strafte 77, 50937 Koln, Germany 
(Dated: February 2, 2008) 

We consider the anisotropic Heisenberg spin- 1/2 chain in a transverse magnetic field at zero 
temperature. We first determine all components of the dynamical structure factor by combining 
exact results with a mean-field approximation recently proposed by Dmitriev et at, JETP 95, 538 
(2002). We then turn to the small anisotropy limit, in which we use field theory methods to obtain 
exact results. We discuss the relevance of our results to Neutron scattering experiments on the ID 
Heisenberg chain compound CS2C0CI4. 
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I. INTRODUCTION 



The spin-i Heisenberg chain is one of the most well- 
understood paradigms for quantum critical behavior. 
The latter term is coined to describe a situation in which 
quantum fluctuations induce critical behavior as a pa- 
rameter like doping or a magnetic field is varied. There 
are many physical realisations of models of (weakly cou- 
pled) Heisenberg chains in anisotropic antiferromagnets. 
Whereas these materials usually display magnetic long- 
range order at very low temperatures T < Tjy, which 
is induced by the weak coupling between chains, there 
is a large "window" above TV where quantum fluctua- 
tions dominate and scaling associated with a quantum 
critical disordered ground state is observed^. An obvi- 
ous question concerns the fate of the intriguing physical 
properties such as the absence of coherent magnon ex- 
citations associated with the "Heisenberg quantum crit- 
ical point", if a magnetic field is applied or "XY"-like 
exchange anisotropics are present. It has been known 
for a long time that critical behaviour persists in both 
cases and there is an entire quantum critical manifold 
rather than just a point. In this regime a host of ex- 
act results on thermodynamic quantities as well as dy- 
namical correlation functions is available (see e.g. Q 
and references therein). On the other hand if one com- 
bines an exchange anisotropy with a magnetic field that 
is at an angle to the anisotropy axis, altogether different 
physics emerges^. The transverse field breaks the contin- 
uous U(l) symmetry of rotations around the anisotropy 
axis and the ground-state develops long-range Neel order. 
Another interesting feature is that the transversality of 
the field prevents the uniform magnetization from reach- 
ing saturation except for infinitely strong fields. 

The main motivation for our work are recent Neutron 
scattering experiments on the quasi one dimensional an- 
tiferromagnet Cs2CoCLf^ in the presence of a magnetic 
field. It was suggested in Ref.Q that an appropriate 



starting point for a description of these experiments is 
the anisotropic spin— 1/2 Heisenberg chain in transverse 
magnetic field 



H = jJ^SJSj+i + &Sj&l+i + SjSj +1 + HSl 



(1) 



with A w 0.25. In view of our interest towards CS2C0CI4, 
we will often consider A = 0.25 in what follows, but our 
results are more general. The inelastic Neutron scatter- 
ing intensity is proportional to 



/(^k)cx^(l 



k a k 



) S*>,k) , 
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where a,0 = x, y, 
S a/3 is defined by 



and the dynamical structure factor 



3,1=1 



dt e-*- R! + ia; * 



x(^(t)Sj_,(0)} 



(3) 



Here the subindex denotes the connected part of the cor- 
relator. For a three dimensional system of uncoupled one 
dimensional chains we have R; = Iciq e, where e is a unit 
vector pointing along the chain direction. The dynamical 
structure factor then depends only on the component of 
the momentum transfer along e, which we will denote by 
k. 

The outline of this paper is as follows: in section ITU wc 
summarize essential steps of the mean-field approxima- 
tion of Dmitriev et al. We then discuss the region of ap- 
plied fields in which we believe the mean-field approxima- 
tion to be applicable. In section lTTll we use exact methods 
to determine all non-zero components of the dynamical 
structure factor within the framework of the mean-field 
approximation. In section Hvl we derive exact results for 
the structure factor in the small anisotropy limit by em- 
ploying field theory methods. In section we discuss a 
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field theory approach to the weak field limit. Section [VI] 
contains a discussion of our results in the context of the 
Neutron scattering experiments on CS2C0CI4. 



II. MEAN-FIELD APPROXIMATION 

Let us first review the main ingredients of the mean- 
field approximation (MFA) proposed by Dmitriev et al. 
in Ref. ja]. We note that the MFA can be applied to 
the XYZ chain in a magnetic field as well. One first 
performs a Jordan- Wigner transformation, which yields 
the spinless fermion Hamiltonian 

n = J2 ^(4^+1 + h - c -) + -f (44+i + h.c.) + 

3 

+ {n 3 - l/2)(n j+1 - 1/2) + H( nj - 1/2). (4) 

Here the hopping and pairing matrix elements are J± = 
(1 ± A)/2 and we have put J — 1. The pairing term is a 
consequence of the anisotropy of the model: it vanishes at 
the isotropic point A = 1. In the MFA the four-fermion 
interaction term is decoupled in all possible ways. The 
relevant expectation values for the on-site magnetization, 
pairing and kinetic term are denoted by 

M = ( nj ) - 1/2, P = (cj +1 Cj), K = (ct +lCi ). (5) 

The mean-field Hamiltonian is then 

V s t ^ — t t 

HmF = 2^ -TT C j C 3 + l + 1T C j C j+l + h - C - 



-1/2), 



(6) 



where we have dropped an unimportant constant, and 
defined the (real) parameters 

J + = J+-2K, J-=J-+2P, H = H + 2M. (7) 

The theory is now purely quadratic in fermion operators, 
and can thus be solved exactly. An alternative picture of 
this theory is obtained by Jordan- Wigner transforming 
back to spin variables, yielding an anisotropic xy model 
in a field, 

^ = E [(1+7)^+1 + (1-7)^+1 + hSj] (8) 
J + 3 

with 7 = J_/J+ and h — H/J+. This model has been 
extensively studied in the literature and many useful re- 
sults on correlation functions are available§iL£i2iiS. For 
our purposes, we first need to complete the mean-field 
setup above by providing the necessary self-consistency 
conditions. First, we use the fact that the free fermionic 
theory can be Fourier transformed and diagonalized using 




FIG. 1: Dispersion relation for A = 0.25 and external mag- 
netic fields H = 0.05, 0.2, 0.4, 0.8, 1.2 and 1.6 (top to bottom 
at k = it. 



a Bogoliubov-de Gennes (BdG) transformation, which 
reads explicitly 



c k = a + (k)c' k + ia_(fc)c'l fe , 
cK = -ia+(k)c' k - a+(k)c'_ k 



with parameters 



a±[k) = T2 



1 ± 



--,1/2 

J+ cos k + H 
uAk) 



(9) 



(10) 



In the thermodynamic limit, the Hamiltonian becomes 



dk 



c'\k) 



u> + (k) 








c'(fc), (11) 



N Jo 2tt~ v ^ 7 V u 
in which the Nambu spinors have bands u>±(k) — ±tu(k), 



uj(k) = J + J (cos k + h) 2 + 7 2 sin k. 



(12) 



This dispersion relation is plotted in Fig for A = 0.25 
and various values of the external field. For A = 0.25 
the gap vanishes for a critical field of H c — 1.604 (cor- 
responding to h = 1). The self-consistency conditions at 
T = can be cast in the form 



dk dui-(k) 
., 2^ dHi 



(13) 



where we have defined the shorthand notations Mq = 
M,Mi = K,M 2 = P,Ho = H,Hi = J+,H 2 = J- 
The numerical solution of these three coupled equations 
determines all the necessary mean-field parameters for 
given values of the anisotropy and of the external mag- 
netic field. 



A. Mean-Field Phase Diagram 

The mean-field phase diagram has already been dis- 
cussed in Ref. . Here we would like to add some obser- 
vations concerning the scaling limits 7 — > and 7 — > 1 
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FIG. 2: Mean-field phase diagram. The critical line is de- 
noted by H C (A) and the classical line by H c i(A). 



to 

U = ~i~ J dx [£ R d x £ R - S, L d x ^ L ] -im J dx £ R £, L , 

(17) 

where £r.l are right and left moving real (Majorana) 
fermions respectively and v — ~ fJ + a ° ; m = -£-(h—l). The 
Majorana fermion theory TFIM is the same as the contin- 
uum limit of the transverse field Ising model (TFIM), i.e. 
the mean- field Hamiltonian for 7=1, with a renormal- 
ized velocity. The elementary excitations are Z2 kinks. 
The region h < 1 corresponds to the ordered phase, 
whereas h > 1 corresponds to the disordered phase of 
the TFIM. In the ordered phase the spin operator a(x) 
has a nonzero expectation value, which corresponds to 
the staggered magnetization in the underlying Heisen- 
berg model. By solving the mean-field equations numer- 
ically we find that 7 1 is only achieved for small A « 
and fields H that are significantly larger than H c . 



respectively. For small 7 and h < 1 we may bosonize 
the mean-field Hamiltonian MFH and obtain a Quantum 
Sine-Gordon model, described by the Hamiltonian den- 
sity 



v 

H =2 



(d x ef + (d x $) 2 } - 7 A cos v^G , (14) 



where $ is a canonical Bose field and its dual field, 
A is a constant and v — J+ v 7 1 — h 2 is the mean-field 
spin velocity at 7 = 0. For h > 1 the 7 = mean-field 
Hamiltonian has a ferromagnetic ground state and an ex- 
citation gap and bosonization cannot be carried out. The 
physics of the Sine-Gordon model SGM is very well un- 
derstood. There are two elementary excitations known 
as soliton and antisoliton. The spectrum consists of scat- 
tering states of an even number of solitons and antisoli- 
tons. By solving the mean-field equations numerically 
we find that small 7 is obtained for small fields H and 
goes to zero as A — > 1 from below. This suggests that 
the physics of the mean-field theory at small fields is ap- 
proximately described by the SGM. We note that for A 
slightly smaller than 1 it was shown in Ref.^lJ that the 
full Hamiltonian Hamil maps onto a Sine-Gordon model. 
We elaborate on this fact in section IIVI For 7 ps 1 it 
is convenient to rewrite the Hamiltonian in terms of the 
real, fermionic operators 



n-l 



n 



/ I I cj)cj' n . 



(15) 



Then one performs a rotation 



t,R,n = 7=—, t,L,n = 7=—, (16) 



V2 



V2 ' 



and finally takes the continuum limit £ Qi „ — > y/ao£ a (x). 
In this way one finds that the Hamiltonian MFH reduces 



B. Staggered Magnetization 

Order in the system is represented by the staggered 
magnetization m s t = ((— lpSJ), which is obtained from 
Ref.fiJ as 



(j + j_) 1 / 4 [i-g 2 /j|] 1/8 

UK, = \ (2(J+ + J_))l/2 





\H\ < J + 
\H\ > h 



(18) 



The mean-field result rhfstag becomes exact on the clas- 
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FIG. 3: Staggered magnetization as a function of external 
field for different anisotropies A = 0.25,0.5,0.75. 
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sical line h c i = \/2(l + A), wherein m st = 
plot m s t as a function of magnetic field for several values 
of A in FigEl 
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C. Applicability of the Mean Field Approximation 

It is immediately clear from the results for the stag- 
gered magnetization shown in Fig|^ that the mean-field 
theory fails to describe the model Hamil in the limit of 
small applied magnetic fields H. In zero field Hamil re- 
duces to a critical XZX spin chain, for which the stag- 
gered magnetization vanishes. The XZX chain has a U(l) 
symmetry corresponding to rotations around the y-axis, 
whereas the MFA breaks this symmetry and is therefore 
invalid. Simple scaling arguments suggest that in the 
limit H — * 0, m s t should be proportional to H a , where 
the exponent a is a function of the anisotropy A 5 . How- 
ever, a renormalisation group analysis of the small H 



J 



regime of Hamil is nontrivial because in a bosonized de- 
scription the transverse field carries conformal spin. We 
review some important features of this analysis in section 

In order to get some rough measure concerning the re- 
gion of magnetic fields H in which the MFA may work 
well, we have computed the one-loop corrections to the 
mean-field values of values Ai, K, P. As there is no 
small parameter in the problem, this calculation should 
be understood as a semiclassical approximation. The cor- 
rections to Ai, K and P are obtained by taking partial 
derivatives of the correction to the ground state energy. 
The latter is found to be 



6En = 



* dpdqdQ 



±F(p+Q,p)F(q+Q,q)cosQ - F(p, q)F(p+Q, q+Q) cos(p-q) F(p+Q,p)F(q+Q, q) cosQ 



* (2tt) 



U)(p) + Lo(q) + w(p + Q)+ uj(q + Q) 
I 



where F{p,q) = ot+(p)et—(q) + a.-(p)a+(q). The correc- 
tions to the mean-field parameters Ai, K, P turn out to 
be quite small except for low fields and near the tran- 
sition: as the gap becomes smaller, fluctuations become 
more important. 

A better way to study the regime of applicability of 
the mean-field approximation is to compare some of its 
predictions to the results of numerical Density Matrix 
Renormalisation Group (DMRG) calculations. Results 
for some values of A have been already been reported by 
Capraro and Gros in Ref. [lflj . 

In Fig0]we compare DMRG results for the magnetisa- 
tion per site with the MFA for A = 0.25. We work with 
periodic boundary conditions and consider system sizes 
of up to 60 sites. The agreement between DMRG and 
the MFA is quite satisfactory, particularly for sufficiently 
large fields H > 1.2J. 

In FigOwe compare DMRG results for the staggered 
magnetisation per site with the MFA for A = 0.25. For 
H < H c the staggered magnetization is computed as 

(o|s*(-ini) , (19) 

where |0) and |1) are the two degenerate ground states 
corresponding to momenta and 7r respectively. It is 
necessary to consider the overlap overlap, because in the 
finite volume and periodic boundary conditions, transla- 
tional invariance implies that 



(o|s:(-irio> = <ii^(-irii> = o 



(20) 



In the infinite volume translational invariance is sponta- 
neously broken and the ground state becomes a linear 
combination of |0) and We note that the finite-size 
effects for weak and strong fields are still rather pro- 
nounced. Such behaviour for the staggered magnetiza- 
tion in finite systems has been previously observed in e.g. 
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FIG. 4: Magnetisation per site calculated from the MFA 
compared to the results of DMRG computations for A = 0.25. 



Ref. . Nevertheless it is clear that the staggered mag- 
netization vanishes above a critical field and approaches 
zero in the weak field limit. The agreement between the 
DMRG results and the MFA for weak fields is as ex- 
pected quite poor. For sufficiently large field H > 1.2J 
the MFA gives good results. Last but not least we have 
used DMRG to determine the gaps of the two lowest 
lying excited states. The results are shown in Fig. 
Both gaps vahish at the critical field H c . For H < H c 
the first excited state is really a degenerate ground state 
with the opposite sign of the staggered magnetization. 
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FIG. 5: Staggered magnetisation calculated from the MFA 
compared to the results of DMRG computations for A = 0.25. 



This degenereacy is removed in the thermodynamic limit 
by spontaneous symmetry breaking. Hence the true ex- 
citation gap is given by the second excited state. We see 



very poor results and fails complately in the limit H — > 0. 
Having in mind these limitations we will now determine 
the dynamical structure factor in the MFA. 



III. DYNAMICAL STRUCTURE FACTOR IN 
THE MEAN-FIELD APPROXIMATION 

The transverse correlators are determined in the frame- 
work of the mean-field theory as follows. Let us denote 
the true ground state by |0) and the mean- field ground 
state by |0)mf- Then we employ the following approxi- 
mation 

(0\SZ(t) S&|0) = (0\e iHat S^e- iHat S&|0) 

-MF^Ie^-'^e-^ S&|0)mf. (21) 

In other words we replace the ground state expectation 
value by the expectation value with respect to the mean- 
field ground state and substitute the xy-Hamiltonian 
MFH for the full Hamiltonian in the time evolution op- 
erator. It is clear that rhfapproxSS will be a good ap- 
proximation as long as the mean-field theory gives an 
accurate description of the ground state and low-lying 
excited states. The transverse spin correlators in the the- 
ory MFH have been determined as a spectral sum over 
intermediate multiparticle states in Refs |?j|3|9|. 
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FIG. 6: Excitation gap calculated from the MFA compared 
to the results of DMRG computations for A = 0.25. 

that the MFA is in very good agreement with the DMRG 
results for fields that are larger than approximately 1.5 J. 
This is only slightly below the critical field H c . 

The above considerations suggest that the MFA works 
very well for both ground state and excited states prop- 
erties as long as the applied field is sufficiently strong 
H > 1.5J. For intermediate field strengths 0.5J < H < 
1.5J we expect the MFA to give at least qualitatively 
correct results. In the small-field regime the MFA gives 



A. Transverse correlations in the "Low-field 
Phase" : h< 1 

For h < 1, the leading contribution comes from two- 
particle intermediate states. As a result all components 
of the structure factor are completely incoherent. The 
other contributions are due to states with 4,6,... par- 
ticles, but they are expected to contribute less spectral 
weight so that we ignore them here. The two-particle 
contributions to the structure factor are given by 



S ap {uo,k) = 

Here 

9xx(ki,k 2 ) - 
g yy (k 1 ,k 2 ) - 

g xv (ki,k 2 ) ■ 



dk\ dk 2 
2tt 



6(lu — u(ki) — uj{k 2 )) x 



xS(k - ki - k 2 + 7r) g a p(ki, k 2 ) . (22) 



c(i + i)- 



[w(fci) - uj(k 2 )Y 



8[1 — cos(fci + k 2 )} uj(ki)uj{k 2 ) 



(7(1 + 7) 1 - cos(fe x - A; 2 ) j 2 



uj{ki)uj(k 2 ) 



+ • 



Cs in(fci^2) 



J. 



u)(ki) oj(k 2 ) 



(23) 



where C = [t?(1 — h 2 )\ 4 and where we have corrected 
a typo in Refill The result for g xy is a conjecture based 
on calculations we have performed in the Ising and Sine- 
Gordon scaling limits and agrees with the known result 



on the classical linei^. We note that the results for 
S aa (uj,k) reduce to the exact results^ 5 for the model 
MFH on the classical line as well, which implies that 
all multiparticle matrix elements vanish. We see that 
the two-particle contribution to the off-diagonal element 
S xy (uj, k) is purely imaginary. This means that it cannot 
be seen in Neutron scattering experiments as it enters 
the expression intensity for the intensity in the form 



(1 



)[s*y(u ) k) + sy*( l j > k)} ) 



(24) 



where we recall that k = k • e is the component of the 
momentum transfer along the chain direction. Using the 
relation 



S yx (uj,k) = {S xy (aj,k)Y 



(25) 



we see that S xy (ui, q) + S yx (u, q) is zero whenever S xy is 
purely imaginary. 

The intensity of the structure factor S aa is plotted 
as a function of k and u> in Figs [7\ to for anisotropy 
A = 0.25 and for external magnetic fields H — 0.8J and 
H = 1.4 J, corresponding to the effective mean- field val- 
ues h = 0.44 and h — 0.82 respectively. It follows from 
our discussion in section III CI that we do not expect the 
results for h — 0.44 to provide a quantitative description 
of the structure factor of the transverse field XXZ chain. 
However the mean-field approximation may still capture 
the redistribution of spectral weight in the various com- 
ponents of the structure factor as the transverse field is 
increased on a qualitative level. One clearly sees in e.g. 
Fig0 that the magnetic field splits the lower boundary 
into two branches. Experimentally, one would thus ex- 
pect to see a double peak in the intensity when scanning 
in frequency for fixed momentum. In Fig. |SJ the evolu- 
tion of the structure factor with increasing magnetic field 
can be seen: the branches recollapse around k = tt, with 
diminishing gap (the gap vanishes at H c = 1.604 J, corre- 
sponding to h = 1). As the field gets closer to the critical 
field, the intensity of the structure factor collapses from 
an incoherent continuum onto an emergent single coher- 
ent mode. 



B. Transverse correlations in the "High-field 
Phase" : h > 1 

For fields h > 1 the dominant contribution to the dy- 
namical structure factor is due to a single-particle coher- 
ent mode with a dispersion relation given by lfP2l . We 
haveZ&a 



S a ^io,k) = f aP {k-ir) 6(u-u(k)) 



(26) 




FIG. 7: Structure factor S xx as a function of k and uj for 
anisotropy A = 0.25 and external magnetic field H — 0.8J. 
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FIG. 8: Structure factor S xx as a function of k and to for 
anisotropy A = 0.25 and external magnetic field H — 1.4J. 
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where 



fxx ^~ [A(k)Y 



fyy(k) ~ Poo ^.(fc)? 
f xz (k) = fyz(k)=0 . 



(27) 



FIG. 9: Structure factor S yy as a function of k and uj for 
anisotropy A = 0.25 and external magnetic field H — 0.8J. 





FIG. 10: Structure factor S vv as a function of k and uj for 
anisotropy A = 0.25 and external magnetic field H — 1.4J. 



FIG. 12: Structure factor S zz as a function of k and u for 
anisotropy A = 0.25 and external magnetic field H — 0.8J. 
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FIG. 11: f xx structure factor amplitude as a function of k 
for external magnetic fields H = 1.604, 1.7 and 2.4 (above the 
critical field). This is the amplitude of the structure factor 
along the dispersion relation 5-function. 




FIG. 13: Structure factor S zz as a function of k and u for 
anisotropy A = 0.25 and external magnetic field H — 1.4J. 



Here we have introduced the parameters 



Al,2 
Poo 

A(k) 



1-7 

[(l-A^i-Ar'xi-A^r 1 )]- 



\\ - 2A 2 cos(fc) + 1 
X[ 2 - 2X[ 1 cos{k) + 1 



(28) 



We note that f xy is purely imaginary, which again means 
that the mixed correlations cannot be observed in Neu- 
tron scattering experiments. 

Near the transition, the amplitude f xx diverges at k = 
7T. For fields higher than the critical field, the divergence 
is smoothed out, as can be seen in figure (|llfl . The gap 
also reopens, as can be seen from equation i|12|) . 



C. Longitudinal correlations 

The longitudinal correlaion function S zz can be com- 
puted directly in terms of a density-density correlator 
of the c fermions. After the BdG transformation to c' 
fcrmions, a straightforward vacuum expectation value 
yields 10 



f(Q + ~,k) = 



k)} 6(w-u(q,k)) 
a(q)a(q + k) - b{q)b{q + k) 



U)(q)uj(q + k) 
tu(q,k) = uj(q- k/2) -Lu(q+k/2)). (29) 

where a(k) — J + cos(fc) + H and b(k) — J_ sin(fc). Equa- 
tion szzmf is valid in both the low-field and in the high- 
field phase. As a result, the longitudinal correlations are 
incoherent for any value of the applied magnetic field. 
The evolution of S zz with increasing applied field H 



8 



is shown in Figs ll2lf3l Firstly, the longitudinal correla- 
tions are generally incommensurate. Second, in low fields 
the intensity is concentrated at the upper boundary of 
the continuum, as is expected as A is small, whereas for 
larger fields a shift towards lower energies is visible. 



IV. FIELD THEORY APPROACH TO THE 
SMALL-ANISOTROPY LIMIT 

The easiest case to deal with by field theory methods is 
the one where the magnetic field H is much stronger than 
the anisotropy 1 — A. The field theory limit in this case 
has been studied in Refill ar± d subsequently in Ref.[l6l|. 
Let us consider the Hamiltonian 

Hzxx,h = J ^ S°~S^ +1 + H ^ Sj 

+ (A~l)Spv +1 =H + H 1 . (30) 

As we assume the field H to be much larger than the 
anisotropy 1 — A, we bosonize at the point A = 1 in the 
presence of a strong field and then switch on the exchange 
anisotropy as a perturbation. 

At low energies Ho is described by a free massless bo- 
son compactified on a ring of radius R i.e. <& and & + 2irR 
are identified. The dual field 9 fulfils 9 = 6 + 1/R. The 
bosonization rules are: 



5/i 



J[x) + (— l) n n(x) + higher harmonics, (31) 



where ag is the lattice spacing, x — nao and 

J x = bcos{f30(x)) sin ~ 26x 

J y = -bsm(/3Q(x)) sin i^jH^) - 2Sx 

J z = Jd x $(x) , 

n x (x) = ccos((3G{x)), n y {x) = csin(/39(x)), 
2ir 

n z (x) = asin(— $(x) - 25x). 
P 



(32) 



Here (3 = 2ttR and the coefficients a, c are known ex- 
actly in the absence of a magnetic field— and numeri- 
cally for several values of the applied field iii2i 9 . The 
(magnetic-field dependent) constant (3, the incommensu- 
rability 5 and spin velocity v s (H) are determined by the 
microscopic parameters A, J and H of the lattice model 
by using the exact Bethe Ansatz solution. This entails 
solving some linear integral equations numerically (see 
e.g. Refs. |l ll20l2 ll22fl ) . The bosonization formulas boso 
hold as long as H < 2 J. For H > 2 J the ground state of 
7io is the saturated ferromagnetic state and the excita- 
tion spectrum is gapped. For later convenience we define 



4(3 2 
8tt - 4/3 2 



(33) 



The perturbing Hamiltonian Tii can now be bosonized 
using boso and fusion of the staggered magnetizations 



gives a contribution proportional to 
cos(2/39) . 



(34) 



In addition there is a small marginal contribution that 
shifts the compactification radius. For simplicity we ne- 
glect it here. 

Thus, at low energies compared to the scale set by the 
applied field H < 2 J, the effective Hamiltonian is given 
by a Sine- Gordon model 



n 



-[(c^) 2 + (a,9) 2 ] - n(A) cos(2/39) . (35) 



The cosine term in the Sine-Gordon model is relevant and 
generates a spectral gap. As 2(3 > ^/An (see e.g. Fig.l 
of Ref . [2fj| ) . the spectrum of the Sine-Gordon model con- 
sists of soliton and antisoliton only. We can immediately 
read off the scaling of the gap as a function of 1 — A 



M oc (1- A)" 



(36) 



The magnetic field dependence enters both via the prefac- 
tor and via the -ff-dependence of (3. In order to calculate 
the prefactor as well as quantities like the magnetization 
we need to know the normalisation of the operator 



(n. — qV qy 



-C cos(2/39) 



(37) 



in the Heisenberg chain in a field, i.e. the Hamiltonian 
hxxz with A = 1. At present the normalisation C is 
not known. In R,ef.|l9l| the issue of how to determine 
C from the large-distance asymptotics of appropriately 
chosen correlation functions in the Heisenberg chain in a 
uniform field has been investigated. In what follows we 
will consider C as a yet unknown function of the magnetic 
field H. The gap is given by 



M 2v T {2^ ) ( (1 - A)Ctt r (1 - 77) \ 2 - 2 " 



2d r (77) 



(38) 



where 

V = /3 2 /2tt , (39) 
and v is the dimensionless spin velocity 

S = fk (40) 

The staggered magnetization in x-direction is nonzero in 
the presence of a transverse field. We have 



((-l)"^)= C (cos/39) , 



(41) 



when 



,23 
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(cos/39) = 



n 

2 




exp 


urn 



smh(rjt) 



2sinh(i) cosh([l - rj\t) 



V - 

r 



(42) 



r 



The magnetization per site can be calculated from the 
ground state energy of the Hamiltonian by taking deriva- 
tives with respect to the magnetic field H. The total 
ground state energy per site is given by 



etot(#) = exxx(ff) + e SG (H, A) , 



(43) 



where the ground state emergy per site of the isotropic 
Heisenberg lattice model exxx(-ff) is known exactly from 
the Bethe Ansatz (see e.g. Refl2lTl and where esc is the 
ground state energy of the Sine-Gordon model24 



e SG (H,A) = - — tan ^ 



(44) 



functions. Let us start with the transverse correlations 
and concentrate on momenta close to n/aQ. The stag- 
gered magnetizations in x and y direction are given by 
boso and calculating their correlation functions reduces 
to the calculation of particular correlators in the SGM. 
This is by now a standard calculation (see e.g. Ref. 20]): 
one determines the first few terms in a Lehmann repre- 
sentation by using the exact matrix elements 2 ^ of the op- 
erator under consideration between the ground state and 
(multi) soliton/antisoliton excited states. The leading 
contribution is due to intermediate states with one soli- 
ton and one antisoliton. Using the notations of Ref.j2(J 
we obtain 



The magnetization is given by 
de tot (H) 



dH 



de SG (H,A) 
dH 



(45) 



The only unknown in the expression for the staggered 
magnetization SGMSMag and magnetization SGMMag 
is the normalization C in normalisation. Once this is 
known with sufficient accurary for taking derivatives with 
respect to H both m and ((— l) n S%) can be evaluated. 



SW{u,A + k)= ^ 

ao ttsVs 2 - 4M 2 



H (-h 



G P G(20 o )/Ci 



£ cosh(^2±HL 



x e 

" \M / 

+contrib. from 4, 6, . . . particles, (46) 
where &h(x) is the Heaviside function, s and 0q are de- 
fined as 



A. Dynamical Structure Factor 



s = ^uj 2 - v 2 k 2 , 6» = arccosh , (47) 



Using the integrability of the SGM it is possible to eval- 
uate the low-energy asymptotics of dynamical correlation 



and Qp and G(9) are given by 



J 



G{9) = iCisinh(0/2)exp 
My/ir r (2=2^) 



Op 



2iJ 



dt /sinh 2 (£[l - iO/ir]) sinh(t[£ - 1]) 
o T V sinh(2t) sinh(£t) cosh(t) 

sinh r\t 



exp 



dt 



t V2sinh(i)cosh([l - rj\t) 2 
I 



V 2 < 



(48) 



where 



The analogous result for S xx (u, ^ + k) is 



C\ = exp 



dtsrah 2 (t/2) sinh(i[£ - 1]) 



t sinh2i sinh£i coshi 



.(49) 



S xx {u : — + k) 
a 



2v,Jc 2 



ttsVs 2 - AM 2 



G(2e )/d 



£ sinh (2^) 



contrib. from 4, 6, . . . particles. 



(50) 
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Finally we determine the longitudinal structure factor 
S zz (uj,k) by making use of recent results by Lukyanov 
and Zamqlodchiko\iSl. It is clear from the bosonization 
formulas boso that the correlations of n z are incommen- 
surate. In other words the longitudinal structure factor 
has low-energy modes at the incommensurate momenta 
± 25. Moreover one can easily establish that the op- 



ao 

erators 



(51) 



have topological charges of =f2 respectively. Here the 
operator of topological charge is defined as 



The leading contribution to the longitudinal spin-spin 
correlation functions is due to intermediate states with 
two solitons or two antisolitons. The formfactors of ver- 
tex operators of the kind VO with definite topological 
charge have been given in Ref.[27|]. A short calculation 
leads to the following result 



S**(w, — ±25 + k) 
ao 



| G (2* )| 2 



4Ws 2 - AM 2 

xe„(-l-2) 



+ contrib. from 4, 6, . . . particles. 



(53) 



Q 







dx d x Q(x). 



(52) The normalization ^(0) is 



Mo) = ^ 



/7r Mr 



JvY 



(i + j) 
(I) 



exp 



dt 

T 



cosh(t) - exp(-[l + 2ir e 



sinh(£t) cosh(i) 



(54) 



where C\ is given in cl. The structure factor depends on 
the transverse field both through the gap M and through 
the parameter 2/3. The variation of 2(3 does not alter 
the various components of the dynamical structure fac- 
tor on a qualitative level: the leading contributions are 
always due to two particles, the structure factor is en- 
tirely incoherent and always vanishes at the threshold 
(as 2/3 is always strictly larger than \/47r), which occurs 
at u> = 2M. However, on a quantitative level the value 
of /3 is quite important: decreasing /3 leads to a narrow- 
ing of the lineshapc in the transverse correlations. What 
precisely me mean by this is illustrated in Fig ll4l where 
we plot J(M/J) 2 -^S vv (uj, ^) as a function of uj/M for 
several values of the magnetic field H . We see that the 
lineshape, when measured in units of the field-dependent 
gap M, sharpens with increasing magnetic field H. How- 
ever, as M itself depends on H, this does not necessarily 
imply that the lineshape measured in physical units like 
meV sharpens with increasing H. 

It is somewhat difficult to compare our field theory 
results to the ones obtained in the mean-field approxi- 
mation as we have no reliable way of determining the 
gap M. 




5 10 15 



co/M 



FIG. 14: JiM/jf-iS^iuJ, ^) as a function of ui/M for H = 
1.373J and H — 1.748J, corresponding to a magnetisation per 
site of 0.2 and 0.3 respectively. The scales of both x and y 
axes depends on H through the gap M. 



V. FIELD THEORY APPROACH TO THE 
WEAK FIELD LIMIT 



The situation where the field is weak compared to the 
anisotropy, i.e. H <g; (1 — A), has been analyzed by renor- 
malisation group methods by Nersesyan et ed in Ref. [2^| . 
It is convenient to rotate the quantization axis such that 



the anisotropy is in z-direction and the field is along the 
x-direction. The effect of the magnetic field is to gen- 
erate an excitation gap and a nonzero expectation value 
for the staggered magnetization in y-direction. This can 
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be seen as follows. The chain Hamiltonian is 

3 

=Ho + H! , (55) 



where Ho is the Hamiltonian of the anisotropic spin-1/2 
chain and 



(56) 



What we do now is to bosonize at the critical point de- 
fined by Ho, i.e. the anisotropic Heisenberg chain, and 
then to perturb away from this fixed point theory by I56|) . 
The bosonized form of Ho is 



Ho 



2 



dx[(d x O) 2 + {d x $) 



(57) 



where $ is a canonical bosonic field and is the dual 
field. The perturbing operator is given by 



Hi = 



bH 

a 



d.r 



cos ( )(-) ) sill ( 



(58) 



where 



(3 = y/2n - 2arccos(A). 



(59) 

The perturbing operator spin is formally relevant, but as 
it is not a Lorentz scalar it requires special treatment. As 
has been shown explicitly for the case A = in Ref.psf 
the Hamiltonian XXZH can be related to a two chain 
model of spinless Luttinger liquids coupled by a weak 
interchain hopping. The generalisation to < A < 1 
is straightforward. In order to match the notations of 
Refl28l we perform a shift 



$(x) 



(60) 



while keeping the dual field unchanged. The 1-loop renor- 
malisation group analysis can then be read off from Refs 
28 29] . At second order in the coupling the following two 
scalar operators are generated radiatively 



c<»s(2 ;<-)). ( ^<J> 



(61) 



Altogether we thus have three perturbing operators with 
corresponding dimcnsionlcss couplings z (of cos(/30) 

cos(^$)), gi (of -cos (jf$)) and g 2 (of - cos(2/39)). 

The RG equations read 28 ' 29 

dz 

m 

d 9i „ f1 / !n 2 

— = 2(1 - -) gi - (77 - -)z 

Cu T) T) 

dg 2 ~.. 1. 



2 j] 



— - 2(1 - 77)52 + ( V - -)z 2 
di J] 



din 77 
dl 



1 



, 2 9i 

= -t:\92V 

2 V r) 



(62) 



where 



(3 2 1 

77(0) = — = 1 arccos(A). 

2tt it 



(63) 



In general these equations have to be solved numerically. 
However, in the case 77(0) w 1, i.e. a small exchange 
anisotropy, one can analyze the equations by a 2-cutoff 
scaling procedure 2 ^. This is because z reaches strong 
coupling while g\ i2 remain small. The gap can then be 
estimated as 



M cx H exp — 



MO) J 



(64) 



The gap is linear in H, but the coefficient is very small 
as 77(0) »1, As is also shown in Ref. 28], at the point Iq 
in the RG flow where z stops remormalizing the Hamil- 
tonian is given by 

h = ^ J dx[(d x e) 2 + (d x $) 2 ] 

3(^871-/77$) - COs(y/8^j&) 



gitv s 

\2iraof 



d.r 



where 



1 - 77" 2 

3 = 2— '— <0 



1 



(65) 



(66) 



As 77 < 1, the first cosine term is irrelevant and can be 
dropped. The 0-field gets pinned at the value which 
implies that the staggered magnetization in j/-direction 



is nonzero 



((-l) n S^) x (sin/36) = const. 



(67) 



If A is not close to 1, the RG equations RG have to 
be integrated numerically. Now several couplings grow 
simultaneously. Eventually 77 becomes very small, while 
gi becomes very large and negative. This means that 
we again flow towards a Sine-Gordon model for the dual 
field and by the same argument as before we find that 
the staggered magnetization in y-direction is nonzero. 



VI. DISCUSSION AND CONCLUSIONS 

As has been discussed recently in Ref.0, CS2C0CI4 is 
a quasi one-dimensional spin-3/2 antiferromagnet with a 
strong single-ion anisotropy D. At energies small com- 
pared to D the spin degrees of freedom are described by 
the anisotropic spin-1/2 Heisenberg model Hamil 4 . Due 
to the smallness of the exchange constant (J = 0.23mcV) 
it is possible to perform inelastic Neutron scattering ex- 
periments in magnetic fields that are fairly large com- 
pared to J. Hence it should be possible to explore much 
of the magnetic phase diagram experimentally. 
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Our analysis suggests that al low fields H < H c , all 
components of the dynamical structure factor are inco- 
herent and the leading contributions come from interme- 
diate states with two particles. These particles are dif- 
ferent from the spinons of the critical XZX chain in that 
they are gapped and do not carry any definite spin (S z is 
not a good quantum number). In a simple picture these 
particles have an interpretation as domain walls between 
the two possible antiferromagnetc spin alignments in the 
x-direction. An argument in favor of this qualitative pic- 
ture may be obtained by consid ering the classical line. 
Here it has been shown in Refs |3ll5j| that the ground 
state is two-fold degenerate. The expectation values of 
the staggered magnetization in the two ground states dif- 
fer by an overall minus sign. This suggests that low-lying 
excitations can be loosely thought of as domain walls be- 
tween the two possible ground states. 

For high fields, the zz-component of the dynamical 
structure factor remains incoherent, whereas the xx, xy 
and yy components now feature a single-particle coher- 
ent mode. For very large fields and small A this particle 
is similar in nature to the Zi kinks of the transverse field 
Ising model in the strong field limit. A physical picture 
for these excitations is presented in Ref.[3(j. 

It would be interesting to compare our results for the 
dynamical structure factor to inelastic Neutron scatter- 
ing experiments on CS2C0CI4 or similar compounds. 

An interesting question which we have not addressed 
is what happens when the magnetic field is applied at 
an arbitrary angle to the exchange anisotropy. The most 
general case is given by 

H = Jj2$j ■ Si+i + (A - 1)S?SV +1 U S. . (68) 
3 



where we may set H x = with out loss of generality (this 
can always be achieved by an appropriate rotation of the 
quantization axis around the y-direction). Rotating the 
spin quantization axis onto the direction of the field leads 
to the following Hamiltonian 

H = JJ2 Sj • S i+ i + H §1 + (A - l)H' jd+1 , (69) 
3 

where H = sgn(H z )|H| and 

H' jd+1 = cos 2 6 S?S? +1 +sm 2 9 S*S] +1 

+ sin 6 cos 6 [5JS? +1 + SjS% +1 ]. (70) 

Due to the presence of the S'JS'J +1 terms in fiprime the 
resulting Hamiltonian can no longer be Jordan- Wigner 
transformed into an expression that is local in terms of 
Fermion operators. Hence the extension of the mean- field 
approximation to the case of an arbitrary orientation of 
the magnetic field is not straightforward. 
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